Dynamic response of artificial bipolar molecules 
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We calculate the equilibrium properties and the dynamic response of two vertically coupled circu- 
lar quantum dots populated by particles of different electrical charge sign, i. e. electrons and holes. 
The equilibrium density profiles are obtained and used to compute the frequencies and oscillator 
strengths of magnetoplasma excitations. We find a strong coupling between the modes derived from 
the center-of-mass modes of the individual dots which leads to an anticrossing with a pronounced 
oscillator strength transfer from the "acoustic" to the "optical" branch. Also, due to breaking of the 
generalized Kohn theorem a number of other than center-of-mass modes are excited whose oscillator 
strengths, however, are rather weak. 

PACS numbers: 73.21. La, 78.67.Hc 
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I. INTRODUCTION 

The physics of quantum dots — the small man-made 
structures in a semiconductor containing anything from a 
few to thousands or more electrons — has already enjoyed 
a busy and fruitful decade.-as.-a major field of research in 
condensed matter physics.BoO A large part of the exper- 
imental work on quantum dots focused on the probing 
of the electronic states insicLe Jig dots by means of far- 
infrared (FIR) spectroscopy.QEfElfl However, a very com- 
mon feature of many quantum structures is their neap] 
circular shape and a parabolic confining potential]' 
which has profound consequences on the optical response. 
According to the generalized Kohn theorem,E3 under 
such conditions the center-of-mass (CM) motion decou- 
ples from the relative motion of the electrons, and the 
electric field of the FIR radiation couples only to the for- 
mer. Consequently, the absorption spectra of circular 
parabolic quantum dots consist only of two CM peaks 
whose positions are independent of the electron num- 
ber and insensitive to the electron-electron interaction 
effects. 

Actual experiments performed on arrays of quan- 
tum dots, containing a fewoO or up to a few hundred 
electrons^ have confirmed the basic two-peak structure in 
the absorption (transmission) spectra. Small deviations 
were explained by taking into account the nonparabolic- 
ity effectSjE^LLateral Coulomb coupling between the neigh- 
boring dots,Eij and spin-orbit interaction.!!] In a separate 
line of development- experimentalists performed a FIR 
spectroscopji-studjO of the formation of incompressible 
edge stripesEij in quantum dots and antidots with a delib- 
erately tailored hard- wall confinement, thereby demon- 
strating that FIR spectroscopy is capable of giving a de- 
tailed insight into many-body systems. As a matter of 
fact, in this work in order to increase the signal strength 
and create non-parabolic potential profiles, double-layer 
dots with three doping layers were used. Therefore, be- 
sides the dominant "optical" modes where the electrons 
in both layers oscillate in phase, also somewhat weaker 
"acoustic" modes were visible. 

Recently, the system of two vertically coupled quantum 



dots was also studied theoretically in greater detail.E3 
This setup is interesting because even in the case 
when the confining potentials of the individual dots are 
parabolic, the Kohn theorem is broken by the interac- 
tion between the non-equivalent dots. Therefore, besides 
the usual CM modes a rather rich spectrum of collective 
modes can be excited, however, most of the oscillator 
strength is still contained in the center-of-mass modes. 

In the present paper we study a similar system con- 
sisting of two vertically coupled circular parabolic quan- 
tum dots containing carriers of opposite charge sign, elec- 
trons and holes, respectively. Our results are very difi 
ferent from those pertaining to the previously studiedta 
vertically coupled dot system and reveal an interesting 
anticrossing of CM-derived modes belonging to sepa- 
rate dots marked by a major oscillator strength trans- 
fer between them. The kind of system we have in mind 
can be structured in bilayer-bipolar heterostructures con- 
taining parallel electron and hole layers in equilibrium. 
These structures have been realized in the crossed gap 
InAs/GaSb systapLZl as well as in biased GaAs/AlGaAs 
heterostructureESEj where electron and hole layers form 
on the opposite sides of an AlGaAs barrier. Most of the 
interest in such systems stems from the possibility (at 
least, in principle) of the formation of Bose-Einstein con- 
densate of indirect excitonsH While the formation of 
the superfluid state has not been demonstrated so far, a 
number of other interesting effects due to the electap.- 
hole coupling have been predicted and/or observedoO 

Our computational approach is based on the general- 
ization to bilayer two-component systems of the formal- 
ism developed by Zaremba and his co-workersE3 andrSjus-. 
cessfully applied to a number of electronic systemsEjca 
This approach is well suited to describe quantum dots 
with a large number of electrons whose dynamic response 
is dominated by collective excitations. lj Our paper has 
the following structure. In Section II we discuss the equi- 
librium charge-density distribution in artificial bipolar 
molecules. The formalism is given in Sec. Ill, and the 



results regarding the dynamic response of these systems 
are presented in Sec. [V. We summarize our results in 
Section [v|. Two Appendices describe, respectively, the 
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calculation of Coulomb integrals and a simplified model 
useful for obtaining quick estimates of the essential char- 
acteristics of the spectrum. 



II. EQUILIBRIUM DENSITIES 

We consider two vertically coupled quantum dots, one 
populated by electrons and the other by an equal num- 
ber of holes. Both electrons and holes are strictly two- 
dimensional (2D) and are laterally confined by parabolic 
potentials. 

The equilibrium and dynamical properties of a many- 
electron system close to the classical regime can be cal- 
culated frpiij— an approximate semiclassical total-energy 
functionaEj'Eil of the electron density n e 



+ / d 2 rw e (r)n e (r) + £ xc [n e (r)]. (1) 



Besides the largest contributions of the direct Coulomb 
interaction energy and the energy in the external confin- 
ing potential w e , the functional (|lj) includes the quantum- 
mechanical kinetic and exchange-correlation energy cor- 
rections. Following previous authors pj'CJ we choose to 
represent the kinetic energy by its lowest-order gradient 
(von Weizsacker) expansion and approximate E xc by the 
local Dirac exchange-only term. We work in the effective 



atomic units defined by setting fi 



e /e = 1. Here 



to* is the effective electron mass (the effective hole mass 
may be different) and e is the dielectric constant of the 
medium. In these units, the kinetic energy and exchange 
functionals, respectively, are given by 



7T f , A f |Vn e (r)| 2 



(2) 



E x [n e ] 
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d 2 r[n e (v)fi\ 



where A = 0.25 is the von Weizsacker coefficient .0 

For a two-component system (electrons and holes) the 
total energy consists of the energies of the two subsystems 
and a coupling term 

E[n ei n h ] = E e [n e ] + E h [n h ] + E coup [n e , n h \. 

The energy functional of the holes Eh[hh] is identical in 
form to that of electrons as given in Eqs. (|l|) and (j|), 
however, since the holes may have a different effective 
mass the kinetic-energy term is scaled by the inverse of 
the ratio of hole-to-electron effective masses k = m* h lm* e . 
In our calculations we use the characteristic value k = 3. 
(Using the bulk GaAs data,c3one obtains k = 7.9 and 1.2 
for heavy and light holes, respectively.) The inter-layer 



coupling is included at the mean-field level 

2 n e (r e )n h (r h ) 



E CO up[n e ,n h ] 



d r e Id r h 



(3) 



with d denoting the vertical separation between the lay- 
ers. 

The equilibrium densities are obtained from the two 
Euler equations 



S 

5n e 
5 
8n h 



E[n e ,n h ] = fi e , 
E[n e ,n h ] = fj, h . 



(4a) 
(4b) 



Following the usual procedureEEEl of evaluating the func- 
tional derivatives in (|J) and expressing them in terms 
of the square root of the particle densities i/) e (^)(r) = 
[^(/^(r)] 1 / 2 we arrive at the equations determining the 
equilibrium charge density profiles 



-■^V 2 + u e (r) - /i £ 



^ e (r) = 0, (5a) 
^(r) = 0. (5b) 
with the following expressions for the effective potentials 



-7T-V 2 +u h (r) - fi h 
2k 



u e (r) 



w e (r) + irip e (r) - \j -il> e {r) 



a r 



- d 



, j>l(r') 



|r-r'- d| 



w h (r) + -7r^(r) - \ -iph(r) 

K V 7T 



2, iPttr') 



d z r 



d 2 r' 



|r-r'- d| 



(6a) 



(6b) 



Since the potentials (|6|) themselves depend on the solu- 
tions of (H), the system of equations (||) and (||) has to be 
solved self-consistently by convergent iterations. The an- 
gular integrations in Coulomb integrals appearing in Eqs. 
(||) and (^) can be carried out analytically as described 
in Appendix The remaining radial equations for i/> e (?i) 
are solved numerically by discretizing the functions and 
potentials on a grid and using an imaginary-time evolu- 
tion technique described in Ref. ^4|. 

The results obtained using the above formalism are 
illustrated in Fig. |l|. In our calculations we use GaAs 
material parameters m* = 0.067 m e and e = 13.4 which 
define the effective Bohr radius (the length unit) a* B = 
% 2 e/m* e e 2 ~ 10 nm and the effective Hartree (the en- 
ergy unit) E* H = e 2 /a* B « 10 meV. Likewise, we express 
frequencies in the units E* H jfi « 1.5 • 10 13 s _1 , and the 
density unit is a* B ^ 2 — 10 16 m~ 2 . In panel (b) of Fig. [J, 
we show two typical examples of the radial distribution of 
the particle densities in the dots obtained by setting the 
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FIG. 1: Equilibrium particle density distributions in bipolar 
quantum-dot molecules. Panel (a) depicts the evolution of the 
particle densities at the center of the dots versus the inter- 
dot distance for equal (full line) and different (dashed line) 
confining potentials. Panel (b) shows the radial dependences 
of the electron and hole densities for the two situations of (a) 
at d = 2. 

confining potentials to uj e — 5, u>h — 3 (dashed line), and 
u} e = ujh = 4 (full line), respectively. The inter-dot sepa- 
ration is d = 2, and each dot contains an equal number 
N e = Nh = 200 of particles. We note that near the center 
the density profiles closely follow the usual semi-elliptic 
shapeja while at the edges the densities are somewhat 
smoothed if compared to the abrupt squaeferoot behav- 
ior predicted by the classical treatment £lli3 This effect 
is due to the included quantum-mechanical corrections 
which make the density to approach-zero asymptotically 
in the classically forbidden region.ES Panel (a) displays 
the particle densities at the central axis of the quantum- 
dot molecule (i.e. the centers of the two dots) as a func- 
tion of the inter-dot distance d. Here we use the same 
values of N e — Nh = 200 and two different sets of con- 
fining frequencies. The full lines are obtained by setting 
uj e and LOh to the same value 4; in this case, due to a 
higher hole effective mass the radius of the dot contain- 
ing holes is smaller and the hole densities at the center 
are considerably larger if compared to the electronic dot. 
The dashed lines illustrate the case of a quantum-dot 
molecule with better balanced radii and densities of its 
two components. Here the confining frequencies are set 
to w e = 5 and u>h = 3. We observe that in all cases the 
particle densities take off rapidly, and hence the dot radii 
shrink, when the inter-dot distance d becomes compara- 
ble or smaller than the dot radii which are typically in 
the range 2 < R < 3.5. 

III. DYNAMIC RESPONSE - THEORY 

When the electron-hole system is perturbed away from 
the equilibrium there develops an internal restoring force. 
The scalar potentials of its components & e (h) acting on 
electron and hole subsystems are given by the functional 
derivatives of the total energy functional with respect to 



the compqnfint densities evaluated at the modified den- 
sity valuesEao 

K(h) = -V* e(Jl ), with * e(h) - E[n e , n h \. (7) 

This internal force, along with the external force due to 
the electric field of the FIR radiation and the Lorentz 
force in the presence of a perpendicular magnetic field, 
enters the set of four linearized hydrodynamic equations: 
the continuity and the force-balance equations for each 
of the two components. However, in order to avoid rep- 
etition, we will write out explicitly and manipulate only 
two generic equations 

^n 1 + V • [n°v] - 0, (8a) 
d 

k— v = -V<f> + ??eE + r/v x w c . (8b) 

The equations describing the hole (electron) layer are ob- 
tained from (||) by setting the charge-sign factor r\ = +1 
(— 1) and supplying an extra subscript h (e) to the den- 
sity, velocity and potential fields. In addition, for elec- 
trons one sets the effective mass factor to k = 1. In 
Eq. (||), the Lorentz force is expressed in terms of the 
vector u) c directed perpendicular to the layers whose ab- 
solute value equals the cyclotron-resonance frequency of 
the electrons. 

The quantities n° and n 1 denote, respectively, the equi- 
librium density and the linear-order deviation. We work 
with stationary fields that depend on time as e~ tuJt and 
consider only dipole excitations of a given circular polar- 
ization. Therefore, in our square-root of density notation 
we express the densities as 

n=(^ + ^ 9 ) 2 , (9) 

here 9 is the angular coordinate. Thus, the two ground- 
state densities equal ^(ft) = ^e(h) an< ^ ^ e nrs t-order 
fluctuations are given by n*, h s = 2ijj e ^ h )(f) e ^e ie < h '> . Ob- 
serve that the fields ip and <j> are both circularly sym- 
metric, and the correct angular dependence of n 1 is ex- 
plicitly included in the factors e ie . The external electric 
field is also taken to be circularly polarized and deriv- 
able from its corresponding scalar potential E = — V<& cx t 
with 3> ext = —Ere lB and E = const. We note that it 
is not necessary to consider also the opposite polariza- 
tion ~ e~ l9 since these results can be obtained from the 
same calculation by simply changing the direction of w c , 
i.e. the sign of its vertical projection. Positive (negative) 
values of oj c correspond to the direction of the electron 
(hole) cyclotron resonance. Substituting the equilibrium 
and oscillating densities from Eq. (||) into Eqs. (||) we 
obtain 

~iLu2^e l6 + div(^> 2 v) = 0, (10a) 
— iuiKV = — V<I> + 77eE + rrw x uj c , (10b) 
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here we also allow for the presence of a small damping 
force, thus making a replacement lu — > Co = lo — ij in Eq. 

fl)- 

The force-balance equation is readily solved by taking 
its cross-product with Co c and us ing the result to eliminate 
the cross-product term in Eq. ( |lOb| ). Straightforward al- 
gebra gives 



/ 2 2 ~ 2\ 

K - k u> )v 



[iKojV® - r/V<I> x Q c ] 
+ [eE x u; c — ir/KxDeE] 



Substituting this expression into the continuity equation 
we will need only the divergence and the radial compo- 
nent of the velocity field. We separate out the angular 
dependence of $ by writing it as $ = f(r)e t6 , and carry- 
ing out the derivative calculations we obtain the wanted 
quantities 

I 1 2 ~ 2\ 
{UJ C — K LO )V r 

= e t9 ^iConf' — r/LUcf- — ir/ConeE + iuj c eE^j , 

and 

(u> 2 - K 2 o) 2 )divv = iCok\7 2 § = iCoKe l6 A r f. 
Here the operator 

d 2 1 d 1 

r dr 2 r dr r 2 

denotes the radial part of the Laplacian. 

Putting everything together we end up with the follow- 
ing generic equation for the charge density fluctuation cj> 



-f~7 - i-r - Cji; [ ftp' + ^A r ftpj -tj c r}f-ip' 



eEtp'(r)U!K — lo c ) 



(11) 



which will generate the two equations for the electron and 
hole layers. We stress that the quantities / (radial parts 
of the potentials $) themselves are linear functionals of 
</>'s. Therefore, the right hand side of Eq. ( pd| ) can be 
represented as a result of an application of a certain linear 
operator on <f>. 

We write the internal potentials <I> e and $h as 



$ e = $ ee + $ e/l , and $ h = <f> h 



hh • 



thus separating the internal potentials created by intcr- 
and intra-layer interactions. The expressions of the re- 
spective contributions are obtained by straightforward 
functional-differentiation from the definition and read 



2 d 



2ir A 

<£hh = —iph^h-Tr^i^hV^h-^hV^h) 

K 2k 



\r - r 

72, 



(12a) 



r — r' 



$eh = -2 / d z r 
§he = -2 I d 2 r 



2 ,^(r')^(r') 



r — r' — d| 



(12b) 
(12c) 
(12d) 



The Coulomb integrals entering the expressions in Eqs. 
( |l2| ) resemble those encountered in the calculation of the 
equilibrium properties in Eqs. (P) and (||). However, in 
the present case we deal with p-wave charge distributions 
of angular dependence ~ e 10 creating p-wave electrostatic 
potentials. The calculation of these integrals is also dis- 
cussed in Appendix |A]. 

For the sake of notational compactness, we introduce 
linear operators C corresponding to the different terms 
in these expressions, so that Eq. (O) can be written as 



-u;(u; 2 c ~Co 2 )^ e +£ ee ^ e +£ eh cj )h = R e , (13a) 
-lo(lo 2 - K 2 Co 2 )(/) h + £ he (f> e + L hh (j) h = R h , (13b) 

where R e = —eEip'(ui + lo c ) and R h = eEt(;'(CoK — lo c ) 
stand for expressions on the right-hand side of Eq. ([ll|) . 

We choose to expand the fields <j) e and 4>h m the set of 
Darwin- Fock functions of the angular momentum M = 1 . 
Thus, 



.AO 



ro 



E 



a n g n {r/r a ), 



ro 



(14a) 
(14b) 



with a n and b n being the expansion coefficients and the 
radial functions 



^(^/^e- 2 / 2ri i (r2) 

are expressed in terms of the associated Laguerre poly- 
nomials L^(x). The expansions ( |l4| ) can be optimized by 
tuning the scaling radius ro, and typically some 20 terms 
are needed in Eqs. ([l4|) to obtain convergent results. 

The remaining task is the numerical calculation of the 
matrix elements of the operators C in the basis ([lj]) lead- 
ing to the coupled set of linear equations 

- lo(lo 2 - Co 2 )a n + ^(£-,a n , + C±,b n ,) = R e n , (15a) 
~lo(lo 2 - n 2 Co 2 )b n + £(££>n' + 4™'M - Rt (15b) 
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FIG. 2: The magnetic field dependence of the four principal 
modes in the spectrum of a bipolar quantum-dot molecule. 
Symbols "+" and "— " indicate the polarisations. Full lines 
correspond to the vertical separation of d — 3, and dashed 
lines depict the decoupled limit d — oo. The insets illustrate 
the relative arrangements of electrical dipoles of the two dots 
pertinent to the respective anticrossing branches. 

which we solve numerically by LU decomposition!! and 
obtain the sets of expansion coefficients a n and b n . This 
enables us to reconstruct the fluctuating charge density 
profiles 4> e (h) in the two layers and evaluate the energy 
dissipation due to the Joule heating 



<2rr ^e(fc)0e(h) 



(16) 



The sum P = P e + Ph determines the absorption rate. 



IV. DYNAMIC RESPONSE 
RESULTS 



NUMERICAL 



Turning to the description of the absorption spectra 
of coupled bipolar quantum dots, we begin by discussing 
the four most conspicuous modes which evolve from the 
centre-of-mass (CM) modes of the two individual dots 
henceforth referred to as CM modes. Later, we proceed 
to describe the higher resonances and low-frequency edge 
modes whose oscillator strengths are inherently weak 
thus rendering them more difficult to observe experimen- 
tally. 

The basic structure of the CM-mode spectrum as a 
function of magnetic field is shown in Fig. |L These 
results are obtained for coupled dots containing N e = 
Nh = 200 particles each with confinement frequencies 



FIG. 3: The magnetic-field dispersion of the anticrossing CM 
modes of negative polarisation plotted for three different val- 
ues of the vertical inter-dot separation d. The symbols "H" 
and "L" label the higher and the lower anticrossing branches. 
Note the widening of the anticrossing gap and its shift towards 
higher frequencies with decreasing d. The inset compares the 
gap-widths obtained from a numerical calculation (triangles) 
and its fit (dashed line) to those obtained from the simplified 
model of Appendix [Em (full line) . 



set to LJ e = 5 and — 3 (as in the second example 
of Sec. II). The vertical separation between the dots 
is d = 3. The encircled symbols "+" and "— " desig- 
nate the directions of the circular polarizations of the 
respective modes. Our convention is to take the direc- 
tion of the electronic cyclotron resonance as "positive" 
("+") and vice versa. The narrow dashed lines indicate 
the positions of the CM modes of decoupled dots (i. e. 
d = oo). In agreement with the generalized Kohn theo- 
rem, at zero magnetic field their frequencies coincide with 
the confinement frequencies oj e and uih, and split into two 
branches at finite magnetic fields. One notes that due to 
the Coulomb coupling between the dots: (i) all the modes 
are slightly displaced upwards with respect to their po- 
sitions at d — oo, and (ii) the two middle modes which 
are polarised in the same "— " direction anticross, while 
the two modes of "+" polarisation reside in distinct fre- 
quency regions and thus interact only very weakly. Both 
points indicate significant differences from the spectra of 
vertically-coupled electronic quantum dotscJ where the 
inter-dot coupling induces shifts of the modes towards 
lower frequencies and no such anticrossing is observed. 
In these systems, the interaction only couples pairs of 
CM modes that both have positive or negative magnetic- 
field dispersion and do not cross in the absence of in- 
teraction. Thus we see that the charge-sign reversal of 
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FIG. 4: The oscillator strengths of the two anticrossing CM 
modes of the negative polarization. The upper (lower) family 
of curves corresponds to the "optical" ("acoustic") branch. 
A pronounced oscillator strength transfer between the modes 
in the magnetic-field range 2 < uj c < 4 is apparent. As in 
Fig. ^, the symbols "H" and "L" indicate the higher and the 
lower anticrossing branches, respectively. The inset depicts 
the results obtained from the simplified model described in 
Appendix Bl The full (dotted) lines correspond to coupled 
(uncoupled) dots. 



particles in one of the dots does indeed induce a sub- 
stantial qualitative difference. In Appendix |B], we show 
that essential features of the CM mode spectrum can be 
captured within a simplified coupled harmonic-oscillator 
model which can be useful in obtaining quick estimates. 

We take a closer look at the anticrossing modes in 
Figs. U and |^ which show, respectively, the behaviour 
of the frequencies of the two anticrossing branches and 
their oscillator strengths for three different values of d. 
We note from Fig. || that as the inter-dot separation be- 
comes smaller and coupling between the dots increases, 
the anticrossing becomes more pronounced while at the 
same time both branches tend to shift towards higher 
frequencies. The inset of Fig. || shows the dependence 
of the size of the anticrossing gap Alo on the inter-dot 
separation d. The triangles denote the actual calculated 
values while the full line shows the result obtained from a 
simple harmonic-oscillator model discussed in Appendix 
[b|. This model predicts that the gap grows proportionally 
to Aui ~ d~ 3 , however, this law is valid only at relatively 
large distances d > 5. We found that for smaller values 
of d the splitting could be reasonably well fitted by the 
dependence Am ~ (d 2 + c?q)~ 3 / 2 indicated by a dashed 
line in the inset of Fig. ||. 

The evolution of the respective oscillator strengths in 
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FIG. 5: FIR absorption (in arbitrary units) in a bipolar 
quantum-dot molecule close to the anticrossing region. We 
show nine absorption curves corresponding to evenly spaced 
magnetic field values and offset vertically by the same amount 
for clarity. Note the "disappearance" of the lower branch. A 
broadening 7 = 0.05 has been used. 



Fig. ^ is rather peculiar and requires a more detailed ex- 
planation. The sum of the oscillator strengths satisfies 
a sum rule, are we normalise it so that their total sum 
equals 1. In the present example the electron and hole 
numbers are set equal while the holes are taken to be 
k = 3 times heavier than the electrons. Therefore, since 
the oscillator strengths scale as ~ N/m, the electronic 
modes possess 3 times higher oscillator strengths. Thus, 
at zero magnetic field the oscillator strengths of the two 
anticrossing modes start from the values close to 0.125 
for the lower-energy mode which is essentially localised 
in the hole subsystem and 0.375 for the higher one. These 
numbers are slightly modified due to the interaction be- 
tween the modes (the role of interaction becomes more 
important at lower values of d) as well as due to the pres- 
ence of other much weaker modes. The rest 50 % of the 
total oscillator strength at oj c — belongs to the modes 
of the opposite "+" polarisation not shown here. 

In the range of magnetic-field strengths 2 < oj c < 4 the 
two modes interact strongly and anticross. In this region, 
one observes a rather pronounced depression in the os- 
cillator strength of the lower branch reaching nearly zero 
value. The missing oscillator strength is transferred to 
the high-frequency branch. This type of behaviour can be 
understood by realising that at the anticrossing point the 
charge density oscillations of the individual dots combine 
together in either "optical" (the two electric dipoles be- 
ing aligned in parallel) or "acoustic" (antiparallel dipoles) 
fashion, as illustrated in the insets of Fig. |^. Naturally, 
the parallel alignment of two dipoles costs more energy, 
and therefore this "optical" mode has a higher oscillation 
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frequency, while at the same time it possesses a larger 
net dipole moment, and consequently, a higher oscilla- 
tor strength. As one notes in the inset of Fig. |], the 
same type of qualitative behaviour is also observed in 
the coupled harmonic-oscillator model of Appendix 
However, as far as oscillator strengths are concerned, its 
quantitative predictions are not trustworthy, and thus 
one can only rely on the more-accurate numerical treat- 
ment. Our calculations predict the above described oscil- 
lations of the oscillator strengths to be quite strong. In 
Fig. U we show a set of absorption lines simuLatiM those 
obtainable in FIR spectroscopy measurementsliltll! 4 ] which 
have been calculated for the case of a vertical separation 
between the dots d = 3. We plot nine lines correspond- 
ing to nine equally spaced values of the magnetic fields 
between lo c = 1 and lj c — 5 (the anticrossing region) 
thus making the "disappearance" of the low-frequency 
branch apparent. Since the fluctuations of the oscillator 
strengths of these modes can be rather large (of order of 
10 % of the total oscillator strength) we expect that the 
above described effect could be readily observed experi- 
mentally. It is worth mentioning that while in the above 
examples we always dealt with equal numbers of elec- 
trons and holes in the dots, the ratio of the electron and 
hole numbers can be useful to balance the distribution 
of oscillator strengths between the modes. For example, 
increasing the number of holes will brighten the modes 
that are mostly due to oscillations in the hole subsys- 
tem thereby compensating for their diminished oscillator 
strengths because of higher hole effective mass. As we 
show in Appendix the strength of the interaction be- 
tween the two dots and the anticrossing gap scale as the 
geometric mean of the electron and hole numbers. 

Besides the above described strong modes the spectra 
of artificial molecules feature a number of other rather 
weak modes that can be classified into higher resonances 
and low-energy edge modes. Due to a stronger locali- 
sation of the charge-density oscillations associated with 
these modes, the number of terms in the expansions Jl4| ) 
has to be increased for an accurate representation. The 
oscillator strengths of these modes, under the conditions 
of our calculations, typically barely reach 10 -4 of the to- 
tal sum of the oscillator strengths. In Fig. ||, we show the 
magnetic field dispersions of the lowest lying high-energy 
modes and several strongest edge modes. These results 
are obtained for the parameter values N e — Nh = 200, 
uo e = 4, uih = 3, and d — 5. The full (dashed) lines de- 
note branches polarised in the positive (negative) direc- 
tion. In general, the locations and dispersions of these 
modes resemble analogous modes as predicted fox. cou- 
pled quantum dots populated solely by electrons.!!!] The 
key difference is that in the present case one can classify 
the modes into those dominated by oscillations in cither 
the electron or hole subsystems. Thus, the two lower 
high-energy resonances in Fig. ^ are mostly due to holes 
while the topmost mode is mainly electronic. The differ- 
ence can be most easily spotted in the relative arrange- 
ment of the positively and negatively polarised branches 




FIG. 6: The magnetic- field dispersion of the higher- 
resonances and edge modes in a bipolar quantum-dot 
molecule. Note that the middle part of the spectrum is cut 
out. The full and dashed lines denote the modes of the posi- 
tive and negative polarisation, respectively. 



at finite magnetic fields. In the electronic mode, the up- 
per branch is polarised in the positive (i. e. electronic 
cyclotron resonance) direction, while in the case of the 
other two high-energy modes the situation is reversed. 
Moreover, the electronic mode features a considerably 
stronger magnetic-field dispersion. This is due to the 
fact that the higher resonances asymptotically approach 
the cyclotron-resonance line which is n = 3 times steeper 
for electrons than for holes. 

The lower part of Fig. ^| shows six most conspicuous 
edge modes of the bipolar quantum-dot molecule. The 
frequency range covered by these modes and their dispe*-, 
sion again resemble the case of electronic coupled dots.t£l 
However, due to the fact that the polarisations of the edge 
modes are determined by the charge-sign of the carriers, 
in the present case the edge modes can be polarised both 
in the cyclotron and the anticyclotron directions. More- 
over, the direction of polarisation of an edge mode also 
betrays the component whose contribution is dominant. 
Thus, the edge modes mostly influenced by the electronic 
subsystem are polarized in the "— " , i.e. anticyclotron di- 
rection of the electrons. These modes a depicted by the 
dashed lines in Fig. 0. Contrariwise, the full lines in Fig. 
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^ denote the modes mostly due to the oscillations in the 
hole subsystem which are polarised in the "+" , i. e. elec- 
tronic cyclotron resonance direction. 

V. SUMMARY 

In conclusion, we made a theoretical investigation of 
the equilibrium density distributions and the far-infrared 
response of bipolar quantum-dot molecules within a hy- 
drodynamic model including the effects due to exchange 
and kinetic- energy in the von Weizsacker approxima- 
tion. The most conspicuous effect we found is the pro- 
nounced anticrossing between two modified center-of- 
mass modes which takes place when the applied magnetic 
field aligns the frequencies of two centre-of-mass modes. 
The additional distinguishing feature of this anticrossing 
is the strongly non-monotonous behaviour of the oscilla- 
tor strengths of the two resulting branches. The oscillator 
strength transfer from the lower "acoustic" to the higher 
"optical" branch is certainly strong enough to be easily 
observable. On the other hand, the higher resonances 
and edge modes which are also excited in the considered 
setup are quite weak. 



Here, t(>o(r) depends only on the radial coordinate r while 
-01 (r) = fi(r)e ie has a p-wave angular dependence, d 
is the vertical separation between the layers which as a 
special case may equal zero. 

The angular integration in of Eq. flAip can be carried 
out analytically in terms of the complete elliptic func- 
tions of the first kind K(k), and the second kind E(k). 
Introducing p 2 = (r + r') 2 + d 2 we obtain 

J^^fdr'r'^lK^y 

(A2) 

J p (r) = ^jfV/M^lc^Y 
with 

C{k) = hE{k)-K{k)]-K{k). 
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APPENDIX A: COULOMB INTEGRALS 



APPENDIX B: COUPLED 
HARMONIC-OSCILLATOR MODEL 



Guided by a similar simple model introduced in Ref. 
PI we show that a number of basic features of the CM- 
mode spectrum can be derived (at least qualitatively) 
from a simplified model featuring two coupled harmonic 
oscillators. The oscillators concentrate the total masses 
and charges of the two coupled dots and interact via the 
potential 



In this problem we encountered two types of Coulomb 
integrals 



J sir) = d 2 r 



J p {r) = 2 Jd 2 r 



r — r' — d| 
Ir-r'-dl ' 



N e N h 



^d 2 + (r e - it,) 2 



d 2d 3 



(r e -r h ) 2 . (Bl) 



Here r e and are the oscillator coordinates, and we used 
(Al) the fact that in our units e 2 /e = 1. Denoting the base 
frequencies of the oscillators by u> e and ujh we write down 
the coupled equations of motion (m* = l,ni£ = k) 



at- , at 2 N e e. N e N h 

N e r e + N e w s e Y e + —— r e xB + — — (r e - r h ) = 0, 



at - .at 2 N h e. N e N h 

N h KT h + N h Kulr h -rjxB (r e - r h ) = 0. 



(B2) 
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Eqs. flB2| ) are solved by introducing the complex variables z e (fo) = x e(h) + iVeVh) an d assuming a harmonic temporal 
dependence 2 e (fc) ~ exp(ia;i). This leads to the following secular equation for the resonance frequencies 



■ U!CU C 
12 , 



N h a 



-uj 2 + u\ — ujuj c /k + N<,VLq/k 



0, 



(B3) 



here we denoted uj c — eB/c and SIq = d~ 3 . From Eq. 
(B3) the CM-mode frequencies can be readily obtained 
as solutions of a quartic equation. However, basing on the 
smallness of the coupling parameter fig, it is possible to 
extract some simpler approximate expressions. We begin 
by noting that the role of the terms proportional to fig 
enterin g th e diagonal and off-diagonal matrix elements 
in Eq. (B3) is diffierent. The former give first-order cor- 
rections to the frequency dispersion, while the latter con- 
tribute to the second order and are important only close 
to coinciding frequencies, thus defining an anticrossing. 

Therefore, we first neglect the off-diagonal perturba- 
tions and solve two decoupled quadratic equations orig- 
inating from the diagonal terms in Eq. (p33|). The solu- 
tions read 



Wl.2 



UJ C 



± 



(B4) 



UJc 

2k 



± 



4k 2 



-N e Ql 



The modes of positive frequ enc ies w\ and 0J3 obtained 
using the upper signs in Eq. (B4) are of "+" circular po- 
larisation while the lower-sign solutions L02 and W4 are 
negative and correspond to modes polarised in the "— " 
direction. We note, that Eq. (B4) predicts that the ab- 
solute values of all four resonance frequencies increase 
due to the inter-dot coupling and gives an estimate of 
the shifts. In the weak coupling regime they grow as 
fig ~ dr z . This conclusion is in agreement with our 
results (see Fig. ||) and underscores a difference of our 
system from electronic double dots considered in Ref. [If]. 
There, all frequency shifts found were negative. 

Assuming that (as in our numerical calculations) to e > 
LOh we find that the modes polarised in "— " direction 
will cross. The crossing point is readily calculated by 
equating L02 = lo^ and equals 



(B5) 



Note, that we present the zero-order solution obtained 
by setting £7q — > which provides an accurate enough 
estimate. Close to this point the off-diagonal perturba- 
tions are important and have to be taken into account 
to introduce an anticrossing behaviour. We use the fact 
that in the vicinity of the resonances of "— " polarisation 
(roots u>2 and L04) the diagonal terms of Eq. (B3) can be 
approximated by 



-(w - <jJx)(oj - LJ2) ' 
-{UJ - U 3 )(UJ - U3i) 



-(uj 2 - wi)(w - L0 2 ), 

-(W4 - W 3 )(w - W4), 



and obtain a quadratic equation 



(UJ - LU 2 )(lJ - U4,)((Jx - UJ2){^3 - W4) 



N e N h 



valid in this region and giving an approximate behaviour 
of the anticrossing modes. This equation can be used, in 
particular, to estimate the size of the anticrossing gap. 
To this end we calculate the difference of its two roots at 
lu c = LUco and L02 = uja, and obtain 



V K(Vl - UJ2){^3 - W4J 



(B6) 



with 0^1,2,3,4 given by Eq. (B4) 



Equations (B4) - (B6) are useful as quick estimates 
of essential features (mode shifts, position and size of 
anticrossing gap) in the CM-mode spectrum of bipolar 
quantum dot molecules. These estimates are obtained 
for weakly coupled dots and can be asked for quantitave 
predictions only in the limit when inter-dot separation 
considerably excee ds t he dot radii. Thus, the gap size 
obtained from Eq. (B6) agrees with the result of accurate 
numerical calculations within 10 % at d = 5. However, 
as we can see in the inset of Fig. || at closer distances the 
gap grows much more slowly than ~ d~ 3 as given by Eq. 
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